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^ Abstract. We present a halo-independent analysis of direct detection data on "light WIMPs," 

y—{ i.e. weakly interacting massive particles with mass close to or below 10 GeV/c^. We include 

^^ new results from silicon CDMS detectors (bounds and excess events), the latest CoGeNT 

■^ acceptances, and recent measurements of low sodium quenching factors in Nal crystals. We 

focus on light WIMPs with spin-independent isospin- conserving and isospin-violating inter- 
actions with nucleons. For these dark matter candidates we find that a low quenching factor 
would make the DAMA modulation incompatible with a reasonable escape velocity for the 
dark matter halo, and that the tension among experimental data tightens in both the isospin- 
conserving and isospin-violating scenarios. We also find that a new although milder tension 
C^ appears between the CoGeNT and DAMA annual modulations on one side and the silicon 

excess events on the other, in that it seems difficult to interpret them as the modulated and 
unmodulated aspects of the same WIMP dark matter signal. 
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1 Introduction 

The nature of dark matter is one of the fundamental unsolved problems of physics and cos- 
mology. Among the dark matter candidates most actively searched for are weakly interacting 
massive particles (WlMPs), which are particles with weakly interacting cross sections and 
masses in the GeV/c^ -10 TeV/c^ range. Of particular interest are WIMPs with mass around 
1-10 GeV/c2, the so-called "light WIMPs." 

Many theoretically-viable light WIMPs with mass < 10 GeV/c^ have been suggested 
since the late 1970s, ranging from heavy massive neutrinos to supersymmetric particles (for 
early papers, see e.g. [1]). Light WIMPs with spin-independent interactions were first singled 
out in 2004 |2] as a viable WIMP interpretation of the DAMA annual modulation [3] com- 
patible with all direct searches at the time. The interest in light WIMPs intensified after the 
first DAMA/LIBRA results of 2008 \4l , and even more in recent years due to other potential 
dark matter hints in the same region. 

Four direct dark matter search experiments (DAMA [5], CoGeNT PH, CRESST-II [9], 
and now CDMS-II [TU]) have data that may be interpreted as signals from light WIMPs. 
DAMA [5] and CoGeNT [^ report annual modulations in their event rates compatible with 
those expected for a dark matter signal [11]. CoGeNT [6l [8], CRESST-II [9j, and very 
recently, CDMS-II [TO] observe an excess of events above their expected backgrounds; the 
excess may be interpreted as due to dark matter WIMPs. In particular, the CDMS-II collab- 
oration, in their latest analysis of their silicon detectors |10i| , has very recently reported three 
events in their WIMP signal region against an expected background of 0.41 events. Data 
were taken between 2007 and 2008 with 8 silicon detectors, in the nuclear recoil energy range 
7 to 100 keV, for a total exposure of 140.2 kg-days prior to the application of selection cuts. 

On the other hand, upper limits have been placed on dark matter WIMPs by other 
direct detection experiments. The most stringent limits on the average (unmodulated) rate 
for hght WIMPs come from the XENONIO [12], XENONIOO [HI [H] and CDMS-II-Ge ex- 
periments [d5j with the addition of SIMPLE [16] in the isospin-violating case. The amplitude 
of an annually modulated signal is directly constrained by CDMS-II-Ge |17j . 

In the following we compare the different experimental results in a way independent of 
the dark halo model |18H20j for dark matter candidates with the standard spin-independent 
WIMP-nucleus interaction, both with isospin- conserving and isospin-violating couplings. 



2 The halo-independent analysis 

We follow the specific method presented in Ref. |20j . which takes into account form factors, 
experimental energy resolutions, acceptances, and efficiencies with arbitrary energy depen- 
dence. 

The differential recoil rate per unit detector mass, typically in units of counts/kg/day/keV, 
for the scattering of WIMPs of mass m off nuclei of mass number A, atomic number Z, and 
mass mA,z is 

dRA,z _ crA,z{E) 
dE 2m/i^ ^ 

where E is the nucleus recoil energy, p is the local WIMP density, p,A,z = fnA,z fn/{mA,z+m) 
is the WIMP-nucleus reduced mass, aA,z{E) is (a multiple of) the WIMP-nucleus differential 
cross-section daA,z/dE = aA,z{E) 'mA,z/'^fJ-\ z'^'^^ ^^^ 

r){v^in,t)= [ ^fc^d^z; (2.2) 

is a velocity integral carrying the only dependence on the (time-dependent) distribution 
/(v, t) of WIMP velocities v relative to the detector. Here 
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is the minimum WIMP speed that can result in a recoil energy E in an elastic scattering 
with the A, Z nucleus. Due to the revolution of the Earth around the Sun, the r] function has 
an annual modulation generally well approximated by the first terms of a harmonic series 

??(l'min, t) = rjoiVrain) + m('Vmin) COs[a;(t - to)] , (2.4) 

where uj = 27r/yr and to is the time of maximum signal. 

For spin-independent interactions (SI), the WIMP-nucleus cross-section can be written 
in terms of the effective WIMP-neutron and WIMP-proton coupling constants /„ and fp as 



<^A,ziE) = <^P^[Z + {A- Z){U/U)? FXz{E) , (2.5) 
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where Up is the WIMP-proton cross-section and F'^^{E) is a nuclear form factor, which we 
take to be a Helm form factor [21j normalized to i^A,z(0) = 1. In most models the couplings 
are isospin-conserving, /„ = fp. Isospin-violating couplings /„ / fp have been considered 
as a possibility to weaken the upper bounds obtained with heavier target elements, which 
being richer in neutrons than lighter elements, have their couplings to WIMPs suppressed 
for/„//p~-0.7[22]. 

It is worth noticing that the function 

r/(Wmin) = -^- Vivmin) (2.6) 

m 

is common to all direct detection experiments. Expressing the data in terms of Wmin and 
^(■ymin), we can therefore compare the different experimental outcomes without any assump- 
tion about the dark halo of our galaxy. More precisely, we can separately compare the 
unmodulated and modulated parts i?o(^min) = o"p(p/m)ryo and 77i(fmin) = crpip/m)r]i. 



Most experiments do not measure the recoil energy E directly, but rather a detected 
energy E' (or a corresponding number of photoelectrons A'pe) subject to measurement un- 
certainties and fluctuations. These are contained in an energy response function Ga,z(-E, E') 
that incorporates the energy resolution aE{E') and the mean value {E') = E Qa^z{E), 
where Qa,z{E) is the quenching factor (the analogous relation for photoelectrons is {Npg) = 
E Q^^(E)). In this context, recoil energies are often quoted in keVnr, while detected ener- 
gies are quoted in keVee (keV electron-equivalent) or directly in photoelectrons. Moreover, 
experiments have an overall counting efficiency or cut acceptance €{E') that depends on E' . 
A compound detector with mass fraction Ca,z in nuclide A, Z has an expected event rate 
equal to 



§j = <E') r dE Y: Ca,z Ga,z{E, E') ^. (2.7) 

•'^ A,Z 

Using dE = (4/i^ zl^^A,z)vmm dvram, we can write the expected rate over a detected energy 
interval [-Ej , E'^ as 
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'E[,E^]= I dVmmT^fi'^^E!^]{Vmm)fi{Vmm)- (2.^ 



Here we have defined the response function for SI WIMP interactions, with 
Ea,z = 2fi\zv'^-^/mA,z, 



-VSI I ^ ^ '^^rmnCA,za%ziEA,z) /"^^ , 

T^iE^^E'^prmn) = }_^ . , ^7^ / dE Ga,z{Ea,z, E ) e{E ). (2. 



Notice that the temporal modulation of Eq. (2.4) produces a modulation of the rate 



R{t) = RQ + Ricos[uj{t-tQ)] (2.10) 

in any energy interval. 

Our task is to gain knowledge on the functions ??o('ymm) and ??i(umin) from measurements 
i?o i Ai?o and Ri it Ai?i of Rq and Ri, respectively, in all the energy intervals probed by 
experiments. This is possible when a range of detected energies [E'i,E'q\ corresponds to 
only one range of i^mm values [fmin.ij Wmin,2]i for example when the measured rate is due to 
interactions with one nuclide only. In this case, [wmin,ii^min,2] is the Wmin interval where the 
response function TZ?^, £;/i('f^min) is significantly different from zero. We approximate this 
interval with fmin,i = Vuiin{E[ - aE{E[)) and 7;min,2 = Vmm{E'2 + cr£;(£'2)) following Ref. [TO]. 
When isotopes of the same element are present, like for Xe or Ge, the i^min intervals of 
the different isotopes almost completely overlap, and we take Wmin,i; ^min,2 to be the Ca,z- 
weighted averages over the isotopes of the element. When there are nuclides belonging to 
very different elements, like Ca and O in CRESST-II, a more complicated procedure should 
be followed (see Ref. ^20j for details). 

Once the [E[,E'2\ range has been mapped to a [timin,!) Wmin,2] range, we can estimate the 
Vmin-weighted averages 
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■^[E[,E'^]— / '^[£;^,£;^](^min) d^min- (2-13) 

3 Experiments considered 

In the following we compare direct detection data in a way similar to Ref. [20j but including 
new data: the new CDMS-II-Si results [10], a new measurement of the Na quenching fac- 
tor [23], and a new measurement of the CoGeNT rate reduction or acceptance factor [8J. In 
this section we summarize the data included in this analysis, stressing the elements which 
are new with respect to Ref. |20j . 

DAM A. We read the modulation amplitudes from Fig. 6 of Ref. |5j. We consider 
scattering off sodium only, since the iodine component is under threshold for low mass WIMPs 
and a reasonable local Galactic escape velocity. We show results for three values of the Na 
quenching factor: QNa = 0.30, Q-^s, = 0.45 (suggested in Ref. pi]), and the energy dependent 
QNa,Coiiar(-E') obtained by interpolating the central values of the data points in Fig. 9 of 
Ref. [23]. No channeling is included, as per Refs. [231125]. 

CoGeNT. We use the list of events, quenching factor, efficiency, exposure times and cos- 
mogenic background given in the 2011 CoGeNT data release [26j. Differently from Ref. [20], 
we use here the acceptance shown in Fig. 20 of Ref. [8], parametrized as C{E) = 1— exp(— ai?), 
with E in keVee and a = 1.32 ("CoGeNT high"), 1.21 ("CoGeNT med."), and 1.10 ("Co- 
GeNT low" ) . The uncertainty of this parametrization is much smaller than the uncertainty 
of the parametrization used in Ref. [20] , which was modeled after Ref. |?7] . As in [5U] , in the 
figures we plot rjQ plus an unknown flat background h^. 

CDMS-II. As in Ref. [20], we use the germanium data (which we call CDMS-II-Ge) 
from the T1Z5 detector for the total rate jl5] . and the results presented in Ref. [17J for 
the modulation analysis. In addition, we also include the recent results from the silicon 
detector analysis in Ref. [lOj, which we call CDMS-II-Si. Since the energy resolution for 
silicon in CDMS-II has not been measured, we use the energy resolution for Ge in Eq. (1) of 
Ref. [28j, a{E) = y^O.2932 + 0.0562£;/keV keV. We corrected a mistake with respect to a{E) 
in Ref. [20j, which however does not change the results of that paper. With three candidate 
events, we calculate the maximum gap upper limit [29j by taking r]{vxain) as a downward 
step function as in Ref. |20| . Assuming the events are a dark matter signal, we bin the recoil 
spectrum in 2 keVnr energy bins, resulting in or 1 events per bin. In the energy bins 
containing one event we use the Poisson central confidence interval of (0.173, 3.30) expected 
events for zero background at the 68% confidence level to draw error bars. 

XENONIOO. We use the last two data releases of Refs. [13] and [H]. The exposure 
in Ref. ^ is 48 kg x 100.9 days, while it is 34 kg x 224.6 days in Ref. [H]. We derive the 
upper limits using the expressions described in Ref. [20j . Due to the two candidate events at 
lower energies in the latter data set, the new limits are less stringent at low energies and more 
stringent at high energies than those derived from the previous data. As in Ref. [20] for the 
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Figure 1. CDMS-II-Si events (crosses), interpreted as a measurement of the unmodulated WIMP 
rate 770, compared to the most stringent upper bounds on 770, for WIMPs with mass m = 9 GeV/c^, 
spin-independent interactions with nuclei, and (left) isospin-conserving fn/ fp = 1 or (right) isospin- 
violating fn/ fp — —0.7 couplings. For the isospin-conserving couplings, there is tension between the 
CDMS-II-Si events and the XENONIOO limits, and a detailed statistical analysis would be necessary 
to quantify their degree of compatibility. For isospin- violating couplings, the CDMS-II-Si events are 
compatible with all bounds (and the most stringent bound between 400 and 800 km/s in iimin is from 
the CDMS-II-Si data themselves). 



energy resolution we use the Poisson fluctuation formula and a Gaussian single photoelectron 
resolution with ctpmt = 0.5 PE fSO]. 

XENONIO. We take the data from Ref. ^^ and compute an upper limit following the 
procedure in Ref. |20| . (We corrected a mistake in the computer program that affected the 
conversion from recoil energy to number of photoelectrons; after this correction, the limit is 
more stringent.) Again for the energy resolution we use the Poisson fluctuation formula and 
single photoelectron Gaussian resolution with (Tpmt = 0.5 PE. The exposure is 1.2 kg xl2.5 
days. We consider the 32 events within the 1.4 keV-10 keV acceptance box in the Phys. Rev. 
Lett, article (not the arXiv preprint). 

SIMPLE. We consider only the Stage 2 p^, a C2CIF5 detector with an exposure of 
6.71 kg days, one observed event above 8 keV, and an expected background of 2.2 it 0.30 
events. We use the Feldman- Cousins method [3l] to place an upper limit of 3.16 expected 
signal events for a 2.2 expected background and 1 observed event, at the 95% confidence 
level. 

CRESST-II. We take the histogram of events in Fig. 11 of Ref. [9j. The electromagnetic 
background is modeled as one e/7 event in the first energy bin of each module. The exposure 
is 730 kg days. We assume a maximum WIMP velocity in the Galaxy such that W recoils 
can be neglected. A light WIMP will scatter however on both Ca and O, thus complicating 
the issue of associating to any range in detected energy [£[,£'2] a range in t;min since the 
correspondence depends on the target mass. To do this, we follow the procedure described 
in Ref. EOl. 



4 Results 

Figs. [1] to [5] show our results for a WIMP with spin- independent interactions and either 
isospin-conserving couplings {fn/ fp = 1) or isospin- violating couplings {fn/ fp = —0.7). In 
all the figures, the vertical axis shows the quantity •fjc'^, which has units of inverse days. 



Despite these units, fjc^ is not the number of WIMPs impinging on the detector per day. 
Notice that rj in the label of the vertical axis stands for either rjo or r]i depending on the 
experiment. In many cases we plot measurements or bounds for both the unmodulated, r]Q, 
and modulated, 771, parts of r/ in the same figure to be able to compare them. In all realistic 
cases one should have rji sufficiently smaller than tjq. 

In Fig. [T]we consider the CDMS-II-Si data points as a measurement of the unmodulated 
rate t/q, and we compare them to the most stringent upper bounds on rjQ for WIMPs with 
mass m = 9 GeV/c^ and spin-independent interactions. The left panel has fn/fp = Ij the 
right panel has fn/fp = —0.7. For the isospin-conserving couplings, there is tension between 
the CDMS-II-Si events and the XENONIO and XENONIOO hmits, and a detailed statistical 
analysis would be necessary to quantify their degree of compatibility. For isospin-violating 
couplings, the CDMS-II-Si events are compatible with all bounds (and the most stringent 
bound between 400 and 800 km/s in Umin is from the CDMS-II-Si data themselves). 

Fig. ^ shows the CoGeNT measurement of (i) the unmodulated part tjq of the velocity 
integral r]{vuim) plus a constant background 60, i.e. rjQ + bQ, and (ii) the modulated part r/i of 
r/(finm)- Both are given as functions of fmin for a WIMP with isospin-conserving couplings 
and mass m = 9 GeV/c^. Three lines are drawn for each r/o + &o cross: one for the central 
value of the acceptance in Fig. 20 of Ref. ^ ("CoGeNT med."), one for the high value of the 
acceptance ("CoGeNT high"), and one for the low value of the acceptance ("CoGeNT low"). 
In comparison with the analogous figure in Ref. [20], Fig. [2] still shows that the modulation 
amplitude of the CoGeNT data (in blue) is large with respect to the average (in brown), 
certainly larger than the few percent modulation in usual halo models. The differences 
between the high, low and medium 770 measurements are now however hardly distinguishable, 
contrary to the situation in Ref. [20|, because the uncertainty in the acceptance factor of 
CoGeNT in Ref. [8] is now smaller than what it was in Ref. [20]. 

In Fig. [3] we compare the measurements of the modulated part of r/(fmin) for m = 9 
GeV/c^ with the CDMS-II-Si data points and the most stringent modulated and unmodulated 
upper bounds. For the modulated part 771 we show the CoGeNT measurements (blue crosses), 
the DAMA measurements (green crosses), and the CDMS-II-Ge bound (magenta line with 
downward arrow). For the unmodulated part r]Q we show the CDMS-II-Si measurements (red 
crosses), the CDMS-II-Ge bound (blue line), the CDMS-II-Si bound (red line), the XENONIO 
bound (orange line), and the XENONIOO bounds (purple lines, dashed for the latest data). 
The three plots on the left column have isospin-conserving couplings fn/fp = 1- The three 
plots on the right column have isospin-violating couplings fn/fp = —0.7. The rows differ 
in the choice of the sodium quenching factor QNa for the DAMA experiment: the first row 
has the large QNa = 0.45 of Ref. [24J, the second row has the medium Qwa = 0.30 used 
by the DAMA/LIBRA collaboration, and the third row has the low and energy-dependent 
QNa,Coiiar(-E') of Ref. [23J. As QNa decreases, the DAMA points (green) move to higher fmin 
(to the right). With the lowest values of <3Na,Coiiar (-£')) close to 0.10 at the lowest energies, 
the higher energy data points of DAMA would be incompatible with reasonable values of the 
Galactic escape velocity, and all the points are rejected by several upper bounds. 

The case of spin-independent isospin-conserving couplings (the left column of Fig. ^ is 
similar to that in Ref. |20j : only the lowest energy bins of the CoGeNT modulation data and 
of the DAMA modulation data with ^Na = 0.45 may escape the XENONIO and XENONIOO 
bounds and the CDMS-II-Ge modulation bound. The situation for the recent CDMS-II-Si 
events appears similar, although the error bars are larger. In this spin-independent isospin- 
conserving case, the unmodulated part 7/0 measured by CDMS-II-Si overlaps the modulated 
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Figure 2. CoGeNT measurements of the quantities 770+^0 and rji as functions of i^min for a WIMP with 
isospin-conserving couphngs /„ = fp and mass to = 9 GeV/c^ . Here bo is a constant background, 770 is 
the unmodulated part of the velocity integral r]{vmin), and 771 is its modulated part. The uncertainty 
in the acceptance, shown here with the three values of 770 "CoGeNT high" , "med." , and "low" , is 
visibly smaller than the uncertainty of the parametrization used in Ref. [5D] . 



part r/i measured by CoGeNT and/or DAMA. The central values of the CDMS-II-Si events 
are in the region excluded by the XENONIOO and CDMS-II-Ge modulation bounds, although 
the error bars of the two lowest energy CDMS-II-Si events extend into the allowed region. 
This makes it difficult to interpret the CoGeNT and DAMA annual modulations on one side 
and the CDMS-II-Si events on the other as the modulated and unmodulated aspects of the 
same WIMP dark matter signal. 

For the case of isospin- violating couplings with fn/ fp — —0.7 (right column of Fig. [s]), 
the CDMS-II-Si upper bounds are now the most restrictive, having taken the role of the 
XENONIOO bounds. This was to be expected, since the particular ratio fn/ fp = —0.7 was 
chosen to decrease the WIMP-Xe couplings. Notice that also the WIMP-Ge and WIMP-Na 
couplings decrease. For example, comparing the left and right columns of Fig. |3j we see that 
the coupHng [Z + {A - Z){fn/ fp)f of the WIMP with Ge in Eq. [2^ decreased, and thus 
77 increased, by a factor of about 500. For Na the factor is smaller, about 70. Thus, with 
fn/ fp = —0.7, the DAMA modulation points are below the CoGeNT data points, while they 
are above for fn/ fp = 1- With isospin-violating couplings the most severe bound come from 
the CDMS-II modulation limit. With respect to the previous analysis in Ref. [20], there 
is now a slight tension between the modulated and unmodulated data also in the case of 
isospin-violating couplings (right column in Fig. pi). The new CDMS-II-Si bounds (red line) 
is lower than the central values of all but the lowest energy modulated CoGeNT data (blue 
crosses), although the CoGeNT error bars extend well into the allowed region. While pushing 
the XENONIO and XENONIOO bounds away from the modulated data, the fn/ fp = -0.7 
isospin-violation pushes at the same time the CDMS-II-Si events below the CoGeNT and 
DAMA modulated data. This makes it difficult, as it was for the isospin-conserving scenario, 
to interpret the CDMS-II-Si events and the modulated direct detection data as coming from 
spin-independent interactions of dark matter WIMPs (unless the CDMS-II-Si events are 
themselves part of a modulated signal) . 

Fig. |4] shows all the measurements and bounds on the modulated and unmodulated 
parts of r]{vja\n) as a function of v-^dn included in this analysis: the CoGeNT and DAMA 
measurements of, and the CDMS-II-Ge bounds on, the modulated part 771; the CDMS-II-Si, 
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Figure 3. Comparison of the modulated and unmodulated measurements of ?7(wj„i„) as a function 
of Winin for TO = 9 GeV/c^. The left column is for isospin-conserving couplings fn/fp = 1; the right 
column is for isospin-violating couplings fn/fp — —0.7. The top ro'w has sodium quenching factor 
in DAMA QNa = 0.45, the central ro-w QNa = 0.30, and the bottom ro-w QNa = QNa,Coi\ai{E) of 
Ref. [23]. The crosses and lines represent: for the modulated part 771, the CoGeNT measurements 
(blue crosses), the DAMA measurements (green crosses), and the CDMS-II-Ge bound (magenta line 
■with downward arrow); for the unmodulated part rjQ, the CDMS-ITSi measurements (red crosses), 
the CDMS-II-Ge bound (blue line), the CDMS-II-Si bound (red hue), the XENONIO bound (orange 
line), and the XENONIOO bounds (purple lines, dashed for the latest data). 



CoGeNT and CRESST-II measurements, and the CDMS-II, XENONIO, XENONIOO and 
SIMPLE bounds on the unmodulated part r/o- In both left and right panels QNa = 0.30 and 
m = 9 GeV/c^. It is assumed that the WIMP couplings are isospin-conserving in the left panel 
an isospin-violating -with fn/fp = —0.7 in the right panel. For isospin-conserving coupling 
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Figure 4. All the measurements and bound on the modulated and unmodulated parts of ??(wmin) as 
a function of Umin: CoGeNT and DAMA measurements of, and CDMS-II-Ge bounds on, the modu- 
lated part r/i; CDMS-II-Si, CoGeNT and CRESST-II measurements, and CDMS-II-Ge, CDMS-H-Si, 
XENONIOO (dashed hue for the latest data), XENONIO and SIMPLE bounds on the unmodulated 
part rio. Here QNa = 0.30, m = 9 GeV/c^, and (left) /„//p = 1, (right) f^/fp = -0.7. In the left 
panel, the XENONIO and XENONIOO bounds exclude all but possibly the lowest energy CoGeNT 
and DAMA bin. For the candidates considered, the CDMS-II-Si measurements (red crosses) either 
overlap (left) or are lower than (right) the CoGeNT and DAMA measurements of the modulated part 
(blue and green crosses, respectively). This makes it difficult to interpret the CDMS-II-Si events as 
the unmodulated part of a WIMP signal, since the modulated part is expected to be smaller than the 
unmodulated part. 
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Figure 5. As in Fig. Wl but for (left) m = 6 GeV/c^ and (right) m — 12 GeV/c^. The main results 
are qualitatively the same as for to = 9 GeV/c^ in Fig. WJ 



the XENONIO bound seems to completely exclude the r]Q + bo measurement of CoGeNT. The 
tension between the two may be alleviated with a smaller rjQ provided it would not conflict 
with the CoGeNT r/i measurement. The CRESST-II measurement of the unmodulated rate 
is incompatible with the CoGeNT and DAMA measurements of the modulated rate for the 
WIMP candidates considered, since the CRESST-II data points either fall above (left) or 
overlap (right) the CoGeNT and DAMA modulation data. 

The last figure, Fig. [5j shows again all the measurements and bounds on the modulated 
and unmodulated parts of ?7(fmin) included in this analysis, but for different WIMP masses 
than in the previous figures, m = Q GeV/c^ in the left panel and m = 12 GeV/c^ in the 



right panel. In both panels we assume isospin-conserving WIMP-nucleon couplings. The 
main results are qualitatively the same as for m = 9 GeV/c^. When m is much smaller than 
the masses of all nuclei involved, the combination mvram becomes independent of m. Thus 
in this limit all data points and bounds presented in the figures scale together in fmin by 
the ratio (9 GeV/mc^). For the target nuclei considered, this argument applies rather well 
to WIMP masses smaller than 9 GeV/c^, with some small deviations for the lightest nuclei 
(O, Na and Si), as a careful examination of Fig. ^ (left) shows. The scaling argument is still 
approximately correct even for somewhat larger masses, such as m = 12 GeV/c^. All the 
data points and limits in Fig. [5] are approximately scaled towards lower or higher fmin values 
with respect to their locations in Fig. Hi those of CRESST-II (oxygen) somewhat less than 
the others. 

5 Conclusions 

We have used the halo-independent method to compare the current direct detection data on 
light dark matter WIMPs (~ 1-10 GeV/c^ in mass). We have followed the specific method 
presented by two of us in Ref. |20j . which takes into account form factors, experimental 
energy resolutions, acceptances, and efficiencies with arbitrary energy dependence. We have 
focused on light WIMPs with spin-independent WIMP-nucleus interactions and with isospin- 
conserving and isospin-violating couplings to protons and neutrons, in the ratios fn/ fp = 1 
and fn/ fp = —0.7, respectively. 

Compared to Ref. [20] , we have added the three events off silicon recently reported by 
the CDMS-II collaboration against an expected background of 0.41 events [TU], together with 
the upper limits implied by the same data. We have also updated the CoGeNT acceptance 
to its latest published form [8] , and have included the case of a very low sodium quenching 
factor QNa for the DAMA experiment as measured in Ref. [23|. As expected, the low QNa 
makes the DAMA modulation data incompatible with a reasonable escape velocity for the 
dark matter halo, and all the points are rejected by several upper bounds. 

For isospin-conserving couplings, similarly to Ref. [20], we conclude that only the lowest 
energy bins of the CoGeNT and DAMA modulation data (the latter with Q^a = 0.45) may 
escape the XENONIO, XENONIOO and the CDMS-II-Ge modulation bounds. For isospin- 
violating couplings fn/ fp = —0.7, a similar but less severe tension has now appeared, because 
the new CDMS-II-Si bounds are below the central values of all but the lowest energy modu- 
lated CoGeNT data, although the CoGeNT error bars extend well into the allowed region. 

The recent CDMS-II-Si events appear in a similar situation of conflict with upper bounds 
if interpreted as due to dark matter WIMPs with spin-independent interactions, although the 
conflict is milder because the error bars are larger. In the isospin-conserving case, the central 
values of the CDMS-II-Si events are in the region excluded by the XENONIOO and CDMS- 
II-Ge modulation bounds, while the error bars of the two lowest energy CDMS-II-Si events 
extend into the allowed region. In the isospin-violating case fn/ fp = —0.7, the CDMS-II-Si 
events correspond to a rate lower than the amplitude of the modulated rates in CoGeNT and 
DAMA. This makes it difficult to interpret the CoGeNT and DAMA annual modulations on 
one side and the CDMS-II-Si events on the other as the modulated and unmodulated aspects 
of the same WIMP dark matter signal. 



10 



Acknowledgments 

P.G. was supported in part by NSF grant PHY-1068111. E.D.N., G.G. and J.-H.H. were 
supported in part by DOE grant DE-FG02-13ER42022. J.-H.H. was also partially supported 
by Spanish Consolider-Ingenio MultiDark (CSD2009-00064). 

References 

[1] B. W. Lee and S. Weinberg, Cosmological Lower Bound on Heavy Neutrino Masses, Phys. Rev. 
Lett. 39 (1977) 165; 

H. Goldberg, Constraint on the Photino Mass from Cosmology, Phys. Rev. Lett. 50, (1983) 
1419 [Erratum-ibid. 103, (2009) 099905]; 

K. Griest, Calculations oj Rates for Direct Detection of Neutralino Dark Matter, Phys. Rev. 
Lett. 61 (1988) 666; 

A. Bottino, F. Donato, N. Fornengo and S. Scopel, Light neutralinos and WIMP direct 
searches, Phys. Rev. D 69 (2004) 037302 [hep-ph/0307303]. 

[2] G. Gclniini and P. Gondolo, DAMA dark matter detection compatible with other searches, 
hep-ph/0405278; 

P. Gondolo and G. Gelniini, Compatibility of DAMA dark matter detection with other searches, 
Phys. Rev. D 71, (2005) 123520 [hep-ph/0504010]. 

R. Bernabei et al. [DAMA CoUaboration] , Dark matter search, Riv. Nuovo Cim. 26N1 (2003) 
1 [astro-ph/0307403]. 

R. Bernabei et al. [DAMA Collaboration], First results from DAMA/LIBRA and the combined 
results with DAMA/Nal, Eur. Phys. J. C 56 (2008) 333 [arXiv:0804.2741 [astro-ph]]. 

R. Bernabei et al. [DAMA and LIBRA Collaborations], A^ew results from DAMA/LIBRA, Eur. 
Phys. J. C 67 (2010) 39 [arXiv:1002.1028 [astro-ph. G A]]. 

C. E. Aalseth et al. [CoGeNT Collaboration], Results from a Search for Light-Mass Dark 
Matter with a P-type Point Contact Germanium Detector, Phys. Rev. Lett. 106 (2011) 131301 
[arXiv: 1002.4703 [astro-ph.CO]]. 

C. E. Aalseth et al. [CoGeNT Collaboration], Search for an Annual Modulation in a P-type 
Point Contact Cermanium Dark Matter Detector, Phys. Rev. Lett. 107 (2011) 141301 
[arXiv:1106.0650 [astro-ph.CO]]. 

C. E. Aalseth et al. [CoGeNT Collaboration], CoCeNT: A Search for Low-Mass Dark Matter 
using p-type Point Contact Cermanium Detectors, arXiv:1208.5737 [astro-ph.CO]. 

G. Angloher et al. [CRESST-II Collaboration] Results from 730 kg days of the CRESST-II 
Dark Matter Search, Eur. Phys. J. C 72 (2012) 1971 [arXiv:1109.0702 [astro-ph.CO]]. 

R. Agnese et al. [CDMS Collaboration] , Dark Matter Search Results Using the Silicon 
Detectors of CDMS II, [arXiv: 1304.4279 [hep-ex]]. 

A.K. Drukier, K. Frcese, and D.N. Spergel, Detecting Cold Dark Matter Candidates, Phys. 
Rev. D D33 (1986) 3495. 

J. Angle et al. [XENONIO Coll.], A search for light dark matter in XENONIO data, Phys. Rev. 
Lett. 107 (2011) 051301 [arXiv: 1104.3088 [astro-ph.CO]]. 



[6- 



\r 



[9 

[lo; 

[11 
[12; 

[13 

[14 



E. Aprile et al. [XENONIOO Collaboration], Dark Matter Results from 100 Live Days of 
XENONIOO Data, Phys. Rev. Lett. 107 (2011) 131302. [arXiv: 1104.2549 [astro-ph.CO]]. 

E. Aprile et al. [XENONIOO Collaboration], Dark Matter Results from 225 Live Days of 
XENONIOO Data, Phys. Rev. Lett. 109 (2012) 181301 [arXiv:1207.5988 [astro-ph.CO]]. 



11 



[15] Z. Ahmed et al. [CDMS-II Collaboration], Results from a Low-Energy Analysis of the CDMS II 
Germanium Data, Phys. Rev. Lett. 106, (2011) 131302 [arXiv: 101 1.2482 [astro-ph.CO]]. 

[16] M. Fclizardo et al. [SIMPLE-II Collaboration], Final Analysis and Results of the Phase II 
SIMPLE Dark Matter Search, Phys. Rev. Lett. 108 (2012) 201302 [arXiv:1106.3014 
[astro-ph.CO]]. 

[17] Z. Ahmed et al. [CDMS Collaboration], Search for annual modulation in low-energy CDMS-II 
data, arXiv:1203.1309 [astro-ph.CO]. 

[18] P. J. Fox, J. Liu and N. Weiner, Integrating Out Astrophysical Uncertainties, Phys. Rev. D 83 
(2011) 103514. [arXiv:1011.1915 [hep-ph]]. 

[19] M. T. Frandscn et al., Resolving astrophysical uncertainties in dark matter direct detection, 
JCAP 1201 (2012) 024 [arXiv: 1111.0292 [hep-ph]]. 

[20] P. Gondolo and G. B. Gelniini, Halo independent comparison of direct dark matter detection 
data, JCAP 1212, (2012) 015 [arXiv:1202.6359 [hep-ph]]. 

[21] R. H. Helm, Inelastic and Elastic Scattering of 187-Mev Electrons from Selected Even-Even 
Nuclei, Phys. Rev. 104 (1956) 1466. 

[22] A. Kurylov and M. Kamionkowski, Generalized analysis of weakly interacting massive particle 
searches, Phys. Rev. D 69 (2004) 063503 [hep-ph/0307185]; 

F. Giuliani, Are direct search experiments sensitive to all spin-independent WIMP candidates?, 
Phys. Rev. Lett. 95 (2005) 101301 [hep-ph/0504157]. 

[23] J. I. Collar, Quenching and Channeling of Nuclear Recoils in NalfTl]: Implications for Dark 
Matter Searches, arXiv: 1302.0796 [physics. ins-det]. 

[24] D. Hooper et al, A Consistent Dark Matter Interpretation For CoGcNT and DAMA/LIBRA, 
Phys. Rev. D 82 (2010) 123509 [arXiv: 1007. 1005 [hep-ph]]. 

[25] N. Bozorgnia, G. B. Gelmini and P. Gondolo, Channeling in direct dark matter detection I: 
channeling fraction m Nal (Tl) crystals, JCAP 1011, (2010) 019. [arXiv:1006.3110 
[astro-ph.CO]]. 

[26] J. Collar, private communication. 

[27] C. Kelso, D. Hooper and M. R. Buckley, Toward A Consistent Picture For CRESST, CoGeNT 
and DAMA, Phys. Rev. D 85 (2012) 043515 [arXiv:1110.5338 [astro-ph.CO]]. 

[28] Z. Ahmed et al. [CDMS Collaboration], Analysis of the low-energy electron-recoil spectrum of 
the CDMS experiment, Phys. Rev. D 81 (2010) 042002 [arXiv:0907.1438 [astro-ph.GA]]. 

[29] S. Yellin, Finding an upper limit in the presence of unknown background, Phys. Rev. D 66, 
(2002) 032005 [physics/0203002]. 

[30] E. Aprile et al. [XENONIOO Collaboration], Likelihood Approach to the First Dark Matter 
Results from XENONIOO, Phys. Rev. D 84 (2011) 052003 [arXiv:1103.0303 [hep-ex]]. 

[31] G. J. Feldman and R. D. Cousins, A Unified approach to the classical statistical analysis of 
small signals, Phys. Rev. D 57 (1998) 3873 [physics/9711021 [physics. data-an]]. 



12 



